Phase interference in antiferromagnetic quantum tunneling with 
an arbitrarily directed magnetic field 

Rong Lii*, Hui Pan, Jia-Lin Zhu, and Bing-Lin Gu 

Center for Advanced Study, Tsinghua University, Beijing 100084, People's Republic of China 

(February 6, 2008) 

Abstract 

The quantum interference effects induced by the topological phase are stud- 
ied analytically in biaxial antiferromagnets with an external magnetic field 
at an arbitrarily angle. This study provides a nontrivial generalization of 
the Kramers degeneracy for equivalent double-well system to coherently spin 
tunneling at ground states as well as low-lying excited states for antiferro- 
magnetic system with asymmetric twin barriers. The spin-phase interference 
effects are found to depend on the orientation of the magnetic field distinctly. 
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One recent experiment on the molecular magnets Fes showed a direct evidence of the 
topological part of the quantum spin phase (Berry phase) in the spin dynamics. The impor- 
tance of the topological Berry phase in spin tunneling was elucidated by Loss et al., and von 
Delft and Henly.i They showed that this term can lead to destructive (for half-integer total 
spins) and instructive (for integer spins) interference between opposite winding tunneling 
paths in single-domain ferromagnetic (FM) particles. While spin-parity effects are some- 
times be related to Kramers degeneracy,!! they typically go beyond the Kramers theorem in 
a rather unexpected way.i Similar spin-parity effects were found in antiferromagnetic (AFM) 
particles, where only the integer excess spins can tunnel but not the half-integer ones.S The 
effects of magnetic field along the hardi and mediumil axis were studied in AFM particles. 

Theoretical result^Hi showed that in AFM particles the exchange energy is enhanced 
to the magnetic anisotropy, which leads to tunneling of the Neel vector a much stronger 
effect than tunneling of magnetization in FM particles. Therefore, the AFM particle is 
expected to be a better candidate for observing quantum tunneling than the FM particle 
with a similar size. Up to now theoretical studies on AFM tunnelin^i^ have been focused 
on spin-phase interference between two opposite ground-state tunneling paths. The spin- 
phase interference between excited-level tunneling paths are unknown for AFM particles. 
Moreover, the previous worksi^ have been confined to the condition that the magnetic field 
be applied along the easy, medium, or hard axis, separately. The purpose of this letter is to 
study the resonant quantum tunneling and spin-phase interference at excited levels in AFM 
particles placed in a magnetic field at an arbitrary angle Ou. Therefore, our study provides 
a nontrivial generalization of the Kramers degeneracy for equivalent double-well system to 
coherently spin tunneling at ground states as well as low-lying excited states for AFM system 
with asymmetric twin barriers caused by the arbitrarily directed magnetic field. 

The system of interest is a small (~ 5nm in radius), single- domain, AFM particle at a 
temperature well below its anisotropy gap. According to the two-sublattice model,S there 
is a strong exchange energy mi ■ 1112/ x± between the two sublattices, where mi and m2 are 
the magnetization vectors of the two sublattices with large, fixed and unequal magnitudes. 



In the following, we assume that mi > m2 and m — rrii — 1712 ^ mi. The system has the 
biaxial symmetry, with x being the easy axis, y being the medium axis, and z being the 
hard axis. The magnetic field is applied in the ZY plane, at an arbitrary angle in the range 
oi < 9h < 7r/2. In order to obtain the tunnel splitting for quantum coherence, we shall 
calculate the path integral: / r'{6'}r'{0} cxp[—Se{0, 0)], where Se is the effective Euclidean 
action for AFM tunneling, 

. . V f j .mi+m2 ( d(f)\ m ( d(f)\ X± ( dOV 

E [9, (j)) — Ea {9, (p) — mHz cos 9 — mHy sin 9 sin 0, and Ea {9, (f)) is the magnetocrystalline 
anisotropy energy V is the volume of the AFM particle, and 7 is the gyromagnetic ratio. 
9 and are the angular components of mi, which can also determine the direction of the 
Neel vector, r = it, and m — hjs/V, where s is the excess spin of the AFM particle due to 
the non-compensation of two sublattices. For this case, the magnetocrystalline anisotropy 
energy is Ea {9, 0) = K± cos^ 9 + K\\ sin^ 9 sin^ 0, where K\\ and K±^ arc the longitudinal and 
the transverse anisotropy coefficients satisfying K± ^ > 0. Therefore, the Neel vector 
is forced to lie in the 9 — 9q plane, and the fluctuations of 9 about ^0 are small. Introducing 
9 = 9q + a, \a\ <^\, E {9, (p) reduces to 

E (a, (p) = K± sin^ 9oa^ + K\\ sin^ 9o (sin — sin 0o)^ 

+2X2 sin 9o cos ^0 (sin — sin 0o)^ oc, (2) 

where cos^o = 'rnHz/2K±, sin0o = mHy/2K\^sin9o — /i sin /y^l — {\h cos 9h)^, A = 
K\\/K^, h = H/Ho, and Hq = 2K\\/m. 

In the semiclassical limit, the dominant contribution to the transition amplitude comes 
from finite action solution (instanton) of the classical equation of motion. Since the config- 
uration space of this problem is a circle, only two types of instantons must be taken into 
account. We use A to denote the instanton passing through the barrier at 9 = 9o, 4> = 7r/2, 

2 



and B as the instanton passing through the barrier at 9 — 9q, (f) — 37r/2. Correspond- 
ingly, there are two kinds of anti-instantons: A~ and B~ . The small barrier at (f) — 7r/2 
is hUs = hU {(j) = 7r/2) = i^||l^sin^ 6*0 (1 — sin0o)^, and the large barrier at = 37r/2 is 
HUl — hU {(f) — 37r/2) = KwYsiri^ Oq (1 + sin0o)^- Performing the Gaussian integration over 
a, we can map the spin system onto a particle moving problem in one-dimensional potential 
well. Now the Euclidean transition amplitude of this system becomes 
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where Stot — 2S — s is the total spins of two sublattices, and S — miV/h'j is the sublattice 
spins. In Eq. (3) we have taken x± = i^\/J-, where J is the exchange interaction between 
two sublattices. The effective mass and effective potential in Eq. (3) are 
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The potential U {(f)) — U {(f> + 2n7r) has an asymmetric twin barrier. The degenerate ground 
states are given by two different types of minima of U {(f)) at 2n7r + 0o and (2n + 1) vr — 0o- 
U {(f)) can be regarded as a superlattice consisting of two sublattices, and the energy spectrum 
can be obtained by applying the Bloch theorem and the tight-binding approximation. The 
translational symmetry is ensured by the possibility of successive 27r extensions. 

The periodic instanton configuration 0^ satisfies the equation of motion ^M. — 
U {(f>p) — —E, where £■ > can be viewed as the classical energy of the pseudoparticle con- 
figuration. Then the periodic instanton A solution for < < C/5 is sin 0^ 
sn(a;ir, ki) is the Jacobian elliptic sine function of modulus ki — / Ti-")(i+/3) 

a = sin0o + ^ KiiFsfn^eo ' ^ = ^^"'^0 " 
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where 



6 = (1 -«)/(! + a), ^1 = ^^o/^i, 



2K\iy /TiM-SuiOq, and (71 = 2/-^(l + a) (1 — {3). The associated Euclidean action is 
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where 



Wa = 4:MUJI 



(5b) 
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Here K (fci), E (ki), and 11 (fci, are the complete elhptic integral of the first, second, and 
third kind, respectively. The similar method can be applied to the periodic instanton B, 
and the resuh is 83 = ^3 + 2El3 for < E < Us, where 



Wa = 4:MUJI 



(6) 



and ^2 = (1 + /3) / (1 - While for Us <E <Ul, the result is 63 = ^3 + 2Ep, where 



W3 = 2M^3 (1 + a) UJ2 



^ Eik,)-^K{k,) + ^^±^l±^n{k,,^,) 



k2 — ^3 ^ ' ^3 ^ ' ^3 



(7) 



with UJ2 = u;o/g2, ^s = {l + P)/{a-P),g2 = p/ (a-/?), and kl = ^^^g^. 

Now we turn to the calculation of level splittings of excited states. For a particle 
moving in a double- well-like potential f/(x), the WKB formula gives the tunnel split- 
tings of the nth degenerate excited levels or the imaginary parts of the nth metastable 
levels as AE„ (or ImK) = exp (-W),0i where (E„) = 27r/t(E„). t (E„) = 



dx 



is the period of the real-time oscillation in the potential well, where 



Xi,2 (En) are the turning points for the particle oscillating inside the potential U (x). For 
the present case, we find that the level splittings for instantons A and B in the domain 



< E < [/5 are A£a(3) = tt:^ exp (-Wa(b)) , where Ia (E) = t3 (E) = {k[), and 



f-AiB) 



k[ = — kf. For Us < E < Ul, the imaginary parts of the metastable energy levels are 
ImE 



exp [-2W3), where (E) 



tK {k'2), and k'2 



1 - kl 



Then we discuss the low energy limit of the level splitting. By using the small oscillator 
approximation for energy near the bottom of the potential well, = (n + 1/2) Q, Q = 



{(PU / dcj)^) ^^^^^ / M. = pK\\y I fiM. sin 6^0 cos 00 , Eqs. (5b) and (6) are expanded as 

Wa(b),„ = WA(ij),o - + ^) + + ^) 
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where "— " for the instanton A, and "+" for the instanton B. Therefore, the tunnel sphttings 
of nth excited levels are found to be 
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Eq. (8b) shows that at finite magnetic field the WKB exponent for instanton A, Wa.o is 
smaller than that for instanton B, Wb,o because the barrier through which instanton B 
must tunnel is higher than that for instanton A. 

It is noted that hASA,n or hASs^n is only the level shift induced by tunneling be- 
tween degenerate excited states through a single barrier. The periodic potential U (0) 
can be regarded as a one-dimensional super lattice with the sublattices A and B. The 
general translation symmetry results in the energy "band" structure which is formally the 
same as that of a one-dimensional tight-binding model in solid state physics, and the en- 
ergy spectrum could be determined by the Bloch theorem. The Bloch states for sublat- 
tices A and B can be written as $a (C; 0) — -^^n^^^'^"VAi(p — (pn), and $b(^,0) = 
^ En e*^('^"+''VB (0 - 0n - a), where 0„ = 2mr + (po, L = N (a + b), a = tt - 20o, and 
b = 7: + 200. Then the total wavefunction (0) is a linear combination of the two Bloch 
states, (0) = qa (0 (C; 0) + <^b (0 (C; 0)- Including the phase contributions of the 
topological term, we derive the secular equation as 



e-i(^-^^)-ASA,n + e'^^-i^)^AEB,n S^-E (0 



aA (0 



0, (10a) 
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and £n = (n + 1/2) The Bloch wave vector ^ = in the first Brillouin zone. Therefore, 
the eigenvalues of Eq. (10a) are 



E^=£n±^J (A^A,n)' + (A^B,n)' + 2 (Af^,^ {/^8B,n) COS 6, (11a) 



where 

e = 27r J s 



1+ U + ^ + ^sin2 0oj cos^o 



+ 25(-V^)coseosin2^o^ (Hb) 



sj \ J 

The tunnel sphtting of nth excited level is 



A£^„ = 2y (A£^A „)^ + (A£^B „)^ + 2 (A£^^^„) (A^B.n) cos 9. We can rederive this result by 
calculating the transition amplitude, or by the effective Hamiltonian method-i 

At zero magnetic field, sin0o = 0, cos 6^0 = 0, A£^A,n = A£^B,n, the tunnel splitting is 
suppressed to zero for the half-integer excess spins, which is in good agreement with the 
Kramers theorem. The presence of a magnetic field perpendicular to the plane of rotation of 
magnetization yields an additional contribution to the topological phase, resulting construc- 
tive and destructive interferences alternatively for both integer and half-integer excess spins. 
Tunneling is thus periodically suppressed. At finite magnetic field, the tunneling spectrum 
of the degenerate nth excited levels depends on the parity of excess spins s. 



A£„ = 2 [{^£A,nf + {^£B,nf ± 2 (A^A.n) (A£b,„) 

, (12) 



X cos 



27r n 1 + ^ + ^sin^^oj s cos 6*0 + 25 (^-^ ) cos6'o sin^ 6^0 



where for integer s, and "— " for half-integer s. The spin-parity effect and the oscillation 
of the tunnel splitting with the filed is shown in Fig. 1. Another important observation is 
that only the z component of the magnetic field (i.e., along the hard axis) can lead to the 
oscillation of the tunnel splitting for the highly anisotropic case. As shown in Fig. 2, even 
a small misalignment of the field with the z axis can completely destroy the oscillation 
effect, and the oscillation is absent when the field is along the medium axis. For small Oh 
the tunnel splitting oscillates with the field, whereas no oscillation is shown up for large 
Oh- In the latter case, a much stronger increase of tunnel splitting with the field is shown. 

6 



This strong dependence on the orientation of the field can be observed for ground-state 
resonance as well as excited-state resonance. As a result, we conclude that the spin-phase 
interference effects depend on the orientation of the external magnetic field distinctly. This 
distinct angular dependence, together with the oscillation of the tunnel splittings with the 
field, may provide an independent experimental test for the spin-phase interference effects 
in AFM particles. In Fig. 3, we plot the tunnel splittings of the ground-state level and 
the first excited level as a function of the magnetic field for integer excess spins. It is 
clearly shown that the splitting is enhanced by quantum tunneling at the excited levels. 
Detailed calculations of the thermodynamic quantities of the system show that the specific 
heats oscillate with the magnetic field and are strongly parity dependent of excess spins at 
sufficiently low temperatures. Due to the topological nature of the Berry phase, these spin- 
parity effects are independent of details such as the magnitude of excess spins, the shape of 
the soliton and the tunneling potential. 

More recently, Wernsdorfer and Sessoli have measured the tunnel splittings in the molecu- 
lar Fcg clusters with the help of an array of micro-SQUIDs, and have found a clear oscillation 
in the tunnel splittings.B Similar spin-phase interference effects observed in ferromagnetic 
Fes cluster are found theoretically in this letter for AFM particles, which may bring a new 
insight to test the manifestation of quantum effects at a macroscopic level and the influence 
of quantum phases on the tunneling bahaviors of spin systems. 
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Figure Captions: 

Fig. 1 The relative tunnel splitting of the first excited level (n = 1) Asi/ AeA,o {H = 0) as 
a function of H/Hq for integer (s = 10) and half-integer (s = 10.5) excess spins at 9h — 0°. 
Here S = 1000, Ki/J = 0.002, and A = K2/K1 = 0.02. 

Fig. 2 The relative tunnel splitting of the first excited level (n = 1) Aei/AeA,o {H = 0) 
as a function of H/Hq for Oh = 0°, 1°, 3°, 5°, 10°, 20° and 90°, respectively. Here s = 10, 
S = 1000, Ki/J = 0.002, and A = K2/K1 = 0.02. 
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Fig. 3 The relative tunnel splitting A£„/ As^^o {H — 0) of the ground-state level {n — 0) 
and the first excited level (n = 1) as a function of H/Hq for integer (s = 10) excess spins at 
9h = 0°. Here S = 1000, KJJ = 0.002, and A = K2/K1 = 0.02. 
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